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PROBLEM 

Extend  the  theory  of  computation  error  generation  in  digital  filters.  Specifically, 
consider  two’s-complement  fixed-point  multiplication  for  digitized  Gaussian  signal  inputs. 


RESULTS 

Deterministic  properties  of  round-off  and  chopping  were  examined  for  two’s- 
complement  fixed-point  multiplication  errors.  It  was  shown  that  a  regularity  exists  in  the 
structure  of  the  errors  as  the  multiplier  input  is  stepped  through  a  series  of  consecutive 
values.  The  error  structure  is  a  function  of  coefficient  value,  number  representation,  and 
word  lengths.  The  error  properties  are  categorized  according  to  families  of  coefficient  values 
denoted  by  the  parameter  v.  This  parameter  is  the  effective  word  length  of  the  error.  This 
regular  structure  allowed  the  use  of  the  Poisson  Summation  Formula  in  deriving  error  statis¬ 
tics  based  on  multiplier  inputs  obtained  from  quantized  random  variables. 

A  finite  impulse  response  filter  structure  was  used  to  define  the  individual  error 
statistics  which  are  used  in  the  calculation  of  the  filter  output  variance  and  spectrum.  The 
error  statistics  which  were  derived  are: 

( 1 )  error  probabilities  (univariate  and  bivariate), 

(2)  mean  error, 

(3)  second  moment  of  the  error, 

(4)  the  cross-correlation  between  a  multiplier  input  and  the  error  for  the  same  or  a 
different  multiplier, 

(5)  the  autocorrelation  of  the  error  for  one  multiplier,  and 

(6)  the  cross-correlation  between  the  errors  for  two  multipliers. 

The  analytical  form  of  each  statistic  was  shown  to  consist  of  two  parts,  an  asymptotic 
part  and  a  decreasing  part.  The  asymptotic  part,  which  is  dependent  on  v,  is  well  behaved 
and  finite  for  finite  filter  input  mean.  The  asymptotic  part  was  shown  to  be  independent  of 
the  quantizer  input  probability  density  function  shape.  The  only  asymptotic  part  which  de¬ 
pends  on  quantizer  input  statistics  appears  in  the  cross-correlation  between  a  multiplier  input 
and  the  error  for  the  same  or  a  different  multiplier.  The  quantizer  input  mean  appears  there. 
The  decreasing  part  is  in  the  form  of  a  summation  which  varies  as  a  function  of  the  input 
standard  deviation.  It  has  the  property  that  it  decreases  to  zero  in  the  limit  as  the  input 
standard  deviation  is  increased.  Specific  forms  of  the  decreasing  part  of  the  statistics  were 
given  for  Gaussian  quantizer  input  probability  density  functions. 


RECOMMENDATIONS 

Examine  the  error  statistics  presented  in  this  report  relative  to  assumptions  and  engi¬ 
neering  guidelines  used  presently.  Show  where  present  guidelines  are  adequate  and  where 
modifications  are  needed. 
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I.  INTRODUCTION 


The  usual  approach  to  the  statistical  analysis  of  the  effects  of  two’s-complement 
(TC)  roundoff  errors  in  digital  filters  makes  use  of  the  following  argument.  Since  multipli¬ 
cation  errors  are  similar  to  analog  signal  quantization  errors,  the  same  statistical  results  are 
assumed  to  hold.  The  multiplication  error  can  be  represented  as  a  white  noise  which  is  zero- 
mean,  which  is  uniformly  distributed  over  the  magnitude  of  the  least  significant  bit  (l.s.b.) 
after  rounding,  and  which  has  zero  cross-correlation  with  the  multiplier  input. 

Chopping  errors  are  also  assumed  uniformly  distributed,  but  they  have  a  non-zero 
mean  for  the  TC  number  representation  [  1  ]-[2) . 

Recent  results  [31  presented  a  more  accurate  model  of  the  error  generation  process 
that  is  valid  for  “large”  standard  deviation  of  the  multiplier  input.  Error  statistics  were 
derived  showing  their  dependence  on  number  representation,  coefficient  value,  chopping  or 
rounding  scheme,  and  certain  statistics  of  the  multiplier  input.  Statistics  for  other  than 
large  standard  deviation  were  not  obtained . 

This  report  derives  analytic  expressions  for  pertinent  statistics  of  TC  multiplication 
errors  for  the  case  of  arbitrary  standard  deviation  of  the  multiplier  input.  It  starts  first  by 
examining  the  form  of  the  finite  impulse  response  (FIR)  filter  and  determining  which  statis¬ 
tics  of  the  error  are  important.  Next,  the  deterministic  TC  error  properties  are  derived.  The 
analog  signal  quantization  process  is  discussed  and  certain  useful  formulas  are  presented 
through  use  of  the  Poisson  Summation  Formula  (hereafter  abbreviated  as  PSF)[4] .  The 
desired  error  statistics  are  then  derived,  also  through  use  of  the  PSF.  Specific  forms  of  the 
statistics  are  presented  for  the  case  of  a  Gaussian  quantizer  input. 


II.  THE  FILTER  MODEL 


This  section  is  presented  as  motivation  for  the  choice  of  error  statistics  which  are 
subsequently  derived  in  this  paper. 

A  block  diagram  representation  of  an  FIR  filter  is  shown  in  Fig.  1  in  the  form  of  a 
tapped  delay  line.  The  number  of  taps  employed  is  equal  to  U  and  the  time  interval  between 
consecutive  data  samples  is  equal  to  T.  The  taps  are  weighted  by  the  filter  coefficient  values 
{a^;  h  =  0, 1,  .  . . ,  U-l } .  These  are  the  quantized  values  which  have  been  determined  through 

some  design  procedure.  The  filter  input  data  are  the  stationary  random  sequence  values 
{  x(iT);  i  =0,±1,±2,  .  .  .}.  This  sequence  is  obtained  by  sampling  and  quantization  (with 
roundoff)  of  the  analog  waveform  x(t).  The  filter  incorporates  the  multiplication  error  as 
the  error  sequence  {e^fi-h]  T);  i  =  0,±1  ,±2,  .  . .  }.  The  index  h  =  0,  1 ,  .  .  . ,  U-l  identifies 

the  tap  the  error  sequence  is  associated  with.  As  shown  in  [3] ,  the  error  sequence  values  are 
deterministically  related  to  the  corresponding  multiplier  input  values.  Hence,  the  time  index 
for  the  error  sequence  is  written  as  shown  to  avoid  confusion  later  on  when  evaluating  the 
error  statistics  associated  with  particular  multiplier  inputs.  The  filter  output  data  are  the 
sequence  values  { y(iT);  i  =  0,  ±1 ,  ±2,  .  .  .  }.  The  number  values  associated  with  x,  a,  and  y 
are  assumed  exactly  representable  by  binary  words  of  lengths  K,  L  and  M  bits,  respectively. 

A  leading  binary  point  is  also  assumed.  These  lengths  are  exclusive  of  sign  bits.  The  input/ 
output  relation  for  this  filter  is 

y(iT)  =  b(iT)  +  w(iT)  (1) 

where 

U-l 

b(iT)=  J  ahx([i-h]T)  (2) 

h=0 

U-l 

w(iT)=  2,  eh([i-h]T)  .  (3) 

h=0 

The  analysis  presented  in  this  paper  does  not  depend  only  on  the  use  of  sampling 
and  quantization  of  an  analog  waveform.  The  filter  input  can  come  from  the  output  of 
another  digital  filter  for  example.  It  turns  out  that  the  statistical  results  can  be  written  in 
terms  of  the  filter  input  alone  or  in  terms  of  the  quantizer  input.  The  latter  is  more  inter¬ 
esting  in  that  specific  forms  for  the  Gaussian  case  can  then  be  applied  to  the  analysis. 

From  [3] ,  multiplier  overflow  occurs  only  for  certain  coefficient  values  and  for  a 
small  number  of  values  of  x,  if  at  all,  at  one  or  the  other  end  point  of  the  range  of  x  (i.e.,  at 
or  close  to  plus  or  minus  1 .0,  depending  on  the  value  of  the  coefficient).  (Values  of  x  and  a 
for  which  overflow  occurs  are  defined  in  this  paper.)  The  sums  represented  by  b  and  w  can 
also  result  in  filter  overflow.  In  the  analysis  it  will  be  assumed  that  the  probability  of  occur¬ 
rence  of  either  kind  of  overflow  is  so  small  as  to  be  negligible.  Practically,  this  can  be  accom¬ 
plished  by  reducing  the  variance  of  the  input  sequence  or  by  increasing  the  input  and  output 
word  lengths  K  and  M. 
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Ihe  sequence  {b(iT)}  in  (2)  is  the  desired  filter  output.  The  sequence  (w(iT)}  repre¬ 
sents  additive  noise  which  is  deterministically  related  to  the  input  sequence.  However,  the 
usual  statistical  descriptions  can  be  employed.  These  are  the  sequence  mean,  the  variance, 
and  the  power  spectral  density. 

The  mean  of  y(iT)  is  given  by  #iy  =  E[y(iT)]  where  E  denotes  expected  value.  It  is 
computed  from  (1)  as 


My  Mx 


U-l  U-l 

i  >h+  • 

h=0  h=0 


The  variance  of  y(iT)  is  given  by  o“  where 

a“  =  E[y2]  (5 

=  ab  +  Ao  ■  (6 

The  desired  variance  is  given  by  og  where 
U-l  U-l 

ol=  J  2  [agahCx([g-h]T)J  .  0 

g=0  h=0 

The  term  Ag  represents  the  change  imposed  on  the  desired  variance  by  the  multiplication 
errors.  It  can  take  on  negative  values  and  is  given  by  the  following  expression: 

U-l  U-l 

A„-  X  S  [2agcVh<[8-»IT)  +  C  ((g^lT)l  .  (S 

g=0  h=0 

The  autocovariance  CvO)  and  cross-covariances  Cv  ,  (•)  and  C,  (•)  are  defined  in  terms 
*  *gfch  g  h 

of  their  corresponding  auto-  and  cross-correlations  and  mean  values  as  follows: 

Cv(.)  =  R*(‘)-Mv>  (5 


cXgeh(*)  RXgeh(,)  “  Mx^eh  * 

CVh(’)sRVh(,)‘#legMeh  ' 


The  power  spectral  density  of  y(iT)  is  given  by  Sy(co)  where 


Sy(co)  =  ^  RyfrjT)  exp(-jcoT?T) 


(12) 


r?=. 


and  Ry(rjT)  is  the  autocorrelation  of  y(iT).  The  autocorrelation  in  turn  can  be  written  as 


Ry(rjT)  =  Rb(rjT)  +  Ar(t?T) 
where  tj  =  0,  ±  1,  ±2 . 


(13) 


The  desired  autocorrelation  is  given  by  Rb(r?T)  where 

U-l  U-l 

RbU?T)  =  2  2  agahRxdg  +  r?~hl T)  •  (14) 

g=0  h=0 

The  term  Ar(tjT)  represents  the  change  in  the  desired  autocorrelation  and,  hence,  the  spec¬ 
trum.  It  is  given  by  the  following  expression: 

U-l  U-l 

Ar(j?T)  =  2  2  [2agRx  <r h([g  +  T?-h)T)  +  R£  c  ((g  +  Tj-hlT)]  .  (15) 

g=0  h=0  8  8 


Thus,  the  spectrum  in  turn  can  be  written  as 
Sy(oo)  =  Sb(cj)  +  As(a>) 


(16) 


where 


Sb(w)  = 


OO 

X  Rb(T?T)  exP(~jwr?T) 

n=~oo 


and 


As(oj)  = 


2 


Ar(tjT)  exp(-jcorjT)  . 


77= -OO 


(17) 


(18) 


The  time  delay  argument  [(g+??-h)T]  comes  about  because  the  same  data  sequence 
x(iT)  is  used  as  input  to  all  multipliers.  Let  xg(iT)  and  xb(iT)  be  the  inputs  to  multipliers 
g  and  h.  However, 
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xg(iT)  =  x([i-g]T) 


xh(iT)  =  x([i  —  h)  T) 


Then,  for  example. 


Rx  =  E[Xg(iT)xh([i +r}]T>] 


=  Rx([g  +  T?-h]T)  .  (21) 

Since  the  multiplication  errors  and  multiplier  outputs  are  deterministically  related  to  the 

inputs,  it  seems  most  natural  to  write  the  second  order  statistics,  Rv  .  and  Rc  ,  ,  in  terms 

gh  egeh 

of  the  time  delays  imposed  on  the  input  sequence  x(iT).  The  approach  to  the  derivation  of 
the  equations  to  be  used  for  these  statistics  depends  on  whether  p  =  1 .0  or  |p|  <  1.0  for  the 
quantizer  input  sequences  used.*  When  g  =  h  in  the  time  delay  argument,  p  =  1 .0  for  p  -  0 
and  |p|  <  1 .0  for  rj  ^  0.  If  g  h,  then  p  =  1 .0  for  h-g  =  17  and  |p|  <  1.0  for  h-g  4=  17.  Note 
that  when  g  =  h  only  one  coefficient  value  is  involved.  Then 

Rx  e,(*)  =  Rx  e  (*?T)  (22) 

xg  h  xgeg 


Re  e  J*)  =  Re  (t?T>  • 

g  h  eg 


When  g  ^  h, 


Rxgeh<->=Rxgeh<(g  +  7?-h]T> 


%eh(,)  =  Regehd8  +  T?-h)T)  •  (25 

Computationally,  (24)  is  no  different  than  (22)  since  only  one  coefficient  value  is  involved. 
However,  (25)  is  not  usually  the  same  as  (23)  since  two  coefficient  values  are  involved.  If 
the  two  coefficients  are  different,  two  different  transformations  from  x  to  e  exist.  Then 
(25)  must  be  used.  If  the  coefficients  are  equal,  (23)  may  be  used.  Note  that  two  separate, 
equal  coefficient  values  lead  to  the  curious  situation  where  the  cross-correlation  of  both 
error  sequences  is  equal  to  the  autocorrelation  of  either  error  sequence  alone  except  for  a 
constant  difference  in  the  time  delay  argument. 


.  is  defined  in  Section  IV. 


p. 


The  statistics  necessary  for  the  evaluation  of  the  covariances  are  then  the  following: 
the  filter  input  mean  jUx,  the  filter  input  autocorrelation  Rx<  * ),  the  mean  multiplication 


error  the  error  autocorrelation  RJ-),  the  error  cross-correlation  R..  e  (•),  and  the 
correlation  between  error  and  multiplier  input  Rx  (■)•  These  are  the  quantities  for 


analytical  expressions  will  be  presented. 


g  h 


cross- 

which 


10 


1 


I 


III.  DETERMINISTIC  ERROR  PROPERTIES 


COMMENTS  ON  NOTATION 

Multiplier  word  lengths  are  defined  as  follows:  input  x,  K  bits  plus  sign;  coefficient  a, 
L  bits  plus  sign;  result  of  roundoff  or  chopping  y,  M  bits  plus  sign;  and  N  bits  which  are 
eliminated  through  roundoff  or  chopping.  The  lengths  K,  L,  M,  and  N  are  arbitrary  constant 
positive  integers  with  the  only  restriction  that  K+L  =  M+N.  The  relationships  among  the 
various  word  lengths  are  shown  in  the  diagram  in  Fig.  2.  The  coefficient  dependent  parame¬ 
ters  v  and  5  are  also  shown. 

The  following  definitions  will  be  used  in  this  report.  Let  k  be  any  positive  or  nega¬ 
tive  integer.  Then  [k]  2N  is  the  set  of  non-negative  integers  such  that  0  <  [k]  -»N  <  (2^-1 ) 
where  [k]  jsj  is  related  to  k  through  the  equivalence  relation 

[kl2N  =  kmod  2N.  (26) 


The  indicator  function  I^tk]  has  the  property  that 


i2nM  " 


1  for  2N-1  <(k]2N<2N-l  , 

. 

0  for  0  <  [k]  ^  <  2N-1  - 1  . 


(27) 


NUMBER  REPRESENTATION 

Consider  the  number  value  b  and  its  (M+N)-bit  binary  fixed-point  machine  represen¬ 
tation  b*.  The  machine  representation  is  treated  as  a  positive  number  and  is  defined  as 

M+N 

b*  =  b0  •  bjb2  •  •  bM+N=  ^  bj2-1  .  (28) 

i=0 

The  following  relationships  exist  between  b  and  b*.  For  non-negative  numbers  b  >  0, 

b*  =  b  .  (29) 

For  negative  numbers  b  <  0, 

b*  =  2  +  b  .  (30) 

The  binary  number  bg  is  called  the  sign  bit  (i.e.,  bQ  =  ( 1  -sgn  b)/2  for  b  ^  0  and  bg  =  0  for 
b  =  0). 
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Let  b  be  written  as 


b  =  u2~M_N  (31) 

where  the  integer  u  can  take  on  the  values 

u  =  0,  1, .  .  .  ,  2M+N-  1  (32) 

for  all  non-negative  numbers  u  >  0,  and  can  take  on  the  values 

U  =  _2M+N  _2  (33) 

for  negative  numbers  u  <  0.  The  remainder  r  is  a  non-negative  number  which  is  defined 
through  the  following.  Let 

M 

b*=  ^  bi2~i  +  r*2_M  (34) 

i=0 


where 


r*  =  r 


~  0  -bM+lbM+2  •  •  •  bM+N 
=  2-N[2M+Nb*]  2N  . 

Let  y*  be  the  machine  representation  of  the  result  of  roundoff  or  chopping  of  b*.  Then 


(35) 


M 


=  2 

i=0 


bj2-1  + 


^M+l^-^  f°r  roundoff, 


0  for  chopping 


(36) 


Note  that  b^j  can  be  written  as 


bM+l  =  I  N12M+Nb*l  ■ 


(37) 


The  value  y  associated  with  y*  can  be  written  as 

y  =  b  +  e  (38) 

where  e  is  the  value  of  the  roundoff  or  chopping  error.  The  relation  between  y  and  y*  is  the 
same  as  that  between  b  and  b*  in  (29)— (30).  Thus,  for  non-negative  numbers  y  >  0, 


y*  =  y 


(39) 


and  for  negative  numbers  y  <  0, 
y*  =  2  +  y  . 


(40) 


ROUNDOFF  AND  CHOPPING 

The  error  can  be  calculated  from  (29) — (38)  as 

e  =  y  -  b 
=  y*  -  b* 


=  _r*2~M  +  < 


b\l+l  2~^  f°r  roundoff, 
0  for  chopping, 


where 


and 


for 


r*2-M  =  2-M-N  [ul  2N 


bM+l2-M  =  2_Ml2N[ul 


u  =  -2M+N . -1,0,  1,...  ,2M+N-1  . 


(41) 


(42) 


(43) 


(44) 


EFFECT  OF  MULTIPLICATION  COEFFICIENT 
Consider  the  product 

b  =  xa 

where 

x  =  k2~K  , 


and 


a  =  C2~l  , 


K  +  L  =  M  +  N  . 


(45) 

(46) 

(47) 

(48) 
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The  integer  k  can  take  on  the  values 

k  =  -2K,  .  .  .  ,  2K-2,  2K-1  .  (49) 

The  intergers  £  are  similarly  described.  (Replace  k  by  8  and  K  by  L  in  (49).)  Then,  from 
(31) 


u  =  k£  . 


(50) 


Consider  values  of  £  of  the  form 

£  =  £'25  (51) 

were  £'  is  a  constant  odd  number  which  can  take  on  the  values 

£' =  ±1,±3,  ...,±(2L_S-1)  (52) 

when  5  =  0,  1,  ...  ,  L-l.  Note  that  £'  can  take  on  the  value  £'  =  -1  when  5  =  L.  The  values 
of  £  of  the  form  (5 1 )  thus  span  the  set  of  all  possible  non-zero  values  of  £  as  defined  by  (47). 
The  values  of  the  error  e  in  (41 )  can  be  evaluated  through  the  substitution 

u  =  k£'25.  (53) 

Simplification  of  the  error  equations  can  be  made  through  introduction  of  the  parameter  v 
where  v  =  N-5.  Thus,  for  example,  for  non-negative  numbers  u  (or  k£'  >  0), 

r*2~M  =  2-M-N[k£'  25)2n 

=  2_M_N25(k£']2N-5  (54) 

=  2-M-"[k£’]2l/ 


and 


bM+12-M  =  2-MI2N[k£'26) 

=  2"M  I2„[k£']  .  (55) 

The  same  results  are  obtained  for  negative  numbers  u  (or  k£'  <  0). 

These  equations  are  valid  when  no  overflow  occurs  and  also  for  the  specified  values 
of  v  and  5.  Overflow  can  occur  for  roundoff  of  positive  numbers  (k£  >  0).  It  cannot  occur 
for  roundoff  of  negative  numbers  (k£  <  0)  for  TC.  Overflow  occurs  when  the  unmodified 
product  b  takes  the  form 
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xa  =  0.1 1  ...  1  lrirj  .  .  .  . 


(56) 


Thus,  the  values  k  that  result  in  overflow  are  those  that  satisfy  the  inequality 

kfi>(2M+,-l)2N-1  .  (5' 

The  effect  of  overflow  is  not  mirrored  in  the  equations  however.  The  interesting  values  of 
v  are  given  by  the  inequalities  1  <  v  <  N  when  L  >  N,  and  (N-L)  <  v  <  N  when  L  <  N. 
Note  that,  when  6  =  L,  the  coefficient  value  is  equal  to  zero  and  no  errors  occur. 


TC  ERROR  PROPERTIES 

In  the  following  sections,  the  integer  form  e  of  the  error  will  be  used.  It  is  defined 
through  (41  H5 1 )  as  the  following: 

e  =  e2M+v 


=  -[kC'l  „  + 


2**  I^IkS']  for  roundoff. 


0  for  chopping  . 


The  following  error  properties  become  evident. 

PROPERTY  1 :  Consider  values  of  k  of  the  form  k  =  k'  +  p2"  where  the  integers  k'  and  p 
are  confined  to  the  ranges 

0  <  k'  <  2"  < 


-2K~v  <  p  <  2K“"  .  (60) 

For  constant  coefficient  a  and  associated  factors  C'  and  v,  all  values  of  k 
with  the  same  constant  k'  map  onto  the  same  value  of  the  error  e.  Further¬ 
more,  the  mapping  from  k'  to  e  is  one-to-one  with  the  mapping  specified 
by  (58). 


relation 


In  much  of  the  following  analysis,  this  property  finds  use  whenever  the  resulting 


e(k')  =  e(k'  +  p21')  (< 

is  employed.  It  is  easily  seen  that  the  range  of  error  values  e  which  can  occur  are  given  as 


(62) 


-2*'-1  +  1  <e  <  2 1,-1 
for  roundoff,  and 

~2V  +1  <  e  <  0  (63) 

for  chopping. 

PROPERTY  2:  Consider  the  two  coefficient  values  a+  >  0  and  a_  <  0  where  a+  =  -a_. 

Also  let  e^(k)  describe  the  mapping  from  k  to  e  for  a^  constant 
where  y  -  +  or  Then  e+(k  >  =  e_(-k)  for  any  integer  k£(  -<*>,  °°).  In 

particular,  e_(k')  =  e+(2l'-k')  for  k'  =  1,2 . 2P-  1 , 

and  e_(k')  =  e+(k')  for  k'  =  0. 

For  constant  coefficient  value,  a  sequence  of  error  values  results  when  the  multiplier 

input  is  stepped  sequentially  through  the  values  k  =  .  .  .  ,  -1 . 0,  1 . This  property  states 

that  the  resulting  sequence  of  error  values  for  a+  is  the  mirror  image  (reflected  around  k  =  0) 
of  the  sequence  that  results  for  a_  =  -a+. 

A  graphical  example  of  the  mappings  from  the  multiplier  input  to  the  error  is  shown 
in  Fig.  3  for  positive  coefficient  values  and  roundoff.  The  patterns  are  grouped  in  columns 
according  to  the  parameter  v  where  v  =  1 ,  2,  3,  4  as  shown.  Only  one  pattern  results  for 
v  ::  I ,  two  for  v  =  2.  four  for  v  =  3,  etc.  Not  all  patterns  for  v  -  4  are  shown.  Other  features 
ot  the  patterns  that  result  are  illustrated  here.  First,  each  possible  error  value  is  mapped  onto 

by  a  set  of  multiplier  input  values  which  differ  by  2V.  Second,  v  describes  the  class  of  coef¬ 
ficient  values  that  result  in  the  same  number  of  possible  error  values  (and  the  same  error 
model  for  large  standard  deviation  of  the  multiplier  input  as  shown  in  (3) ).  Third,  the 
mapping  from  the  multiplier  input  to  the  error  value  is  dependent  on  the  value  of  the  coeffi¬ 
cient  for  a  given  v.  And  fourth,  for  some  word  length  combinations,  the  same  mapping  can 
occur  for  two  different  coefficient  values  for  a  given  v. 

The  patterns  that  occur  for  chopping  are  the  same  as  for  roundoff  for  each  coefficient 
value.  The  differences  lie  in  a  shift  both  in  the  horizontal  (multiplier  input  variable)  and 
vertical  (multiplication  error)  direction.  Thus,  the  same  comments  hold  regarding  the  v 
class,  number  of  error  values,  reflection  and  other  properties. 
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IV.  MULTIPLIER  INPUT  STATISTICS 


In  this  section  the  analog  signal  quantization  process  is  reviewed.  Relevant  quantizer 
output  (hence  filter  input)  statistics  are  presented  ending  with  the  Gaussian  case.  Statistics 
are  the  filter  input  mean  and  the  input  auto-  and  cross-correlation  for  zero  and  non-zero 
time  lag.  Use  of  the  PSF  is  demonstrated  and  the  notation  applied  to  the  rest  of  the  paper  is 
discussed. 


UNIVARIATE  CHARACTERISTICS 

Let  x  be  a  discrete  random  variable  (r.v.)  which  is  obtained  from  a  continuously  dis¬ 
tributed  r.v.  x  by  quantization  with  roundoff.  That  is,  let  x  =  kq  where  the  integer  k  is 

^  1/ 

chosen  such  that  q(k-Vi)  <  x  <  q(k+‘/2)  and  q  =  2  ^  is  the  quantization  step  size.  The  A/D 
converter  used  has  an  output  word  length  of  K  bits  plus  sign  with  TC  number  representation. 
This  restricts  the  values  of  the  A/D  converter  output  to  correspond  to  the  range  denoted  by 

k2~^  where  k  =  -2^,  -2^  +1 , ....  2^-1 .  In  the  following,  it  is  mathematically  convenient 
to  assume  x  (and  hence  k)  is  not  bounded.  In  order  to  practically  realize  this  assumption, 
the  A/D  converter  input  can  be  made  small  enough  in  most  cases  so  that  the  effect  of  quan¬ 
tizer  overflow  on  the  derived  statistics  is  negligible. 

In  the  following  it  is  assumed  that  x  is  statistically  stationary  and  has  a  probability 
density  function  (p.d.f.)  fj(x)  associated  with  it.  This  p.d.f.  is  also  assumed  to  be  continu¬ 
ous  everywhere  on  the  real  line.  The  probabilities  of  occurrence  of  each  k  are  written  as 
Px(k)  and  are  defined  as 

r  q(k+!/2) 

Px(k)=  /  f~(x)dx  (64) 

J  q(k-'/2) 

for  k  =  0,  ±  1 ,  ±2 . The  dependence  on  K  is  assumed  understood. 

Associated  with  the  p.d.f.  f^(x)  is  the  characteristic  function  Q^fcu)  defined  by  the 

Fourier  transform 


Q*<w)  = 


fx(£)  exp(jw£)d£  . 


(65) 


The  special  case  that  will  be  used  is  that  of  the  Gaussian  p.d.f.  f^(x)  where 


with  characteristic  function  given  by 


Q?(«)  =  exp 


2  2 
o  OJ 


w  w  , 

- 2 —  +JUOJ 


(67) 


The  expression  Px(k)  is  well-behaved  for  any  real  k.  Its  Fourier  transform  Qx(w)  is 


given  by 


°x(W)=/- 


Pv(£)exp(jw£)d{ 


(68) 


=  sinc(o)/2ir)Q^(co/q) 

where  sinc(x)  =  sin(jrx)/irx. 

Note  that  [4] 

(&QHu=<rWq  <69) 

and 

Rjj(0)  -  •  (70) 

cj=0 


MEAN  MULTIPLIER  INPUT 

The  mean  multiplier  input  nx  can  be  obtained  as  follows: 

Mx  =  E[x] 

=  2~KE[k]  (71) 

OO 

=  q  2  kPx(k) . 

k=-oo 

Since  the  function  kPx(k)  is  well  behaved  for  any  real  k,  the  PSF  can  be  employed  on  (71 ) 
to  yield 
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A 


V 


OO 


S=— oo 


cj=2jts 


=  ?x  + 


where,  from  (69), 


and 


S  Q;(2ws/q)  • 


s#0 


(72) 


(73) 


(74) 


Note  that,  aside  from  stationarity  and  continuity  assumptions,  no  special  form  of  f~(x)  is 

assumed  at  this  point.  The  Gaussian  form  of  i*ix  is  obtained  through  substitution  of  (67) 
into  (74).  The  result  is 


OO 

W  =  S  ^T^exp  {“-(ws)2(a/q)2 J  sin(27rsp/q)  . 
s=l 


(75) 


NOTATION  AND  ASSUMPTIONS 

There  are  some  features  of  ~jlx  and  |/ix  which  are  present  in  subsequent  equations 
where  the  arrow  notation  is  used.  As  can  be  seen  from  (75)  the  function  fjux  is  dependent 

on  two  parameters,  p  and  o,  each  in  relation  to  the  quantization  step  size  q.  The  function 
can  be  made  to  approach  arbitrarily  close  to  zero  through  a  suitable  increase  in  the  parame¬ 
ter  a  for  any  constant  value  of  ft. 

Mathematical  limits  of  #ix  can  be  taken  in  either  one  of  two  ways.  The  first  way  is 
to  increase  a  and  keep  n  constant.  The  second  way  is  to  increase  a  and  keep  the  ratio  o/n 
constant.  The  first  limit  has  the  distinction  that  /ix  is  a  constant.  For  the  second  limit  7rx 

increases  as  o.  It  will  be  assumed  that  the  horizontal  arrow  denotes  a  variable  which  is  either 
constant  or  depends  in  a  well-behaved  manner  on  n  and/or  a  for  a  finite  n  and  a.  The 
vertical  arrow  denotes  a  variable  which  tends  to  zero  as  a  increases  regardless  of  the  assump¬ 
tions  about  fi. 

As  with  the  mean  quantizer  output,  subsequent  cases  will  occur  where  a  progression 
is  made  from  a  statistic  expressed  as  an  expected  value  to  the  asymptotic  (->•)  and  decreasing 
(i)  terms.  It  will  be  assumed  without  comment  that  the  functions  involved  are  well-behaved 
and  that  the  PSF  was  used  to  arrive  at  the  final  results. 
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MULTIPLIER  INPUT  SECOND  MOMENT 


The  second  moment  of  the  multiplier  input  is  Rx(0).  It  is  the  autocorrelation  func¬ 
tion  for  zero  time  lag.  It  can  be  obtained  as  follows: 

R„(0)  =  E[x21 


=  2~2K  E[k2] 


=  q2  ^  k2Px(k) 

k=-oo 


(76) 


=  -q2 


I 


s=- 


du> 


2Qx(<j) 


=  tx(0)  +  ;Rx(0) 


6J=27TS 


where,  from  (70) 

Sx(0)=R3f(0)  +  q2/12 

=  a2  +n2  +  q2/ 1 2  , 
and,  for  the  Gaussian  case 


(77) 


OO 

=  ^  (-1 


|Rx(0)  =  2q^  Z.  (“I  f 
s=l 


-Us  +  2(a/q)2 


2(7TSr 


'  |  ^  (n/ q)  ]  sin(27rs(n/q)) 


cos(27rs#x/q) 


exp(-2(7rs)2(a/q)2) 


(78) 


BIVARIATE  CHARACTERISTICS 

By  analogy  with  (64),  the  joint  probability  of  occurrence  of  a  pair  of  multiplier 
inputs  (k  j  ,k2)  is  written  as  Px  ^  x^(k  j  ,k2)  which  can  be  specified  in  the  same  way  as  the  set 

of  probabilities  (Px(k)}.  Let  the  multiplier  input  r.v.’s  x  j  and  x2  be  obtained  from  the 
joint  r.v.’s  Xj  and  x2  by  quantization  with  roundoff.  Assume  x  j  and  x2  are  stationary  and 
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have  a  joint'p.d.f.  that  is  continuous  everywhere  in  the  interval 

-°o<  xj  ,x->  <  Further  assume  x  j  and  \2  are  correlated  with  normalized  correlation 

coefficient  p  =  ^(^/(o^)  wkere  Ipl  <  1  and  O]  =  and  a2  =  o y  .  Thus 

r  q(k,+4)  /*q(k 2+^) 


/•^lT  x. 

PxjX2(kl’k2J=  /  .  / 


f  J  q(k2  - 


,  f  x ,  x2C^l  2)  (Jx  1  dx  2 


for  ki  ,k-)  =  0,  ±1 ,  ±2 . 


Associated  with  the  p.d.f.  pXj *s  t^ie  characteristic  function  j  (upa^) 
defined  by  the  Fourier  transform 


Qxj^l’^)3 


Ilf- 


fxjX2^1’^2^exP^£jJl^l  +jw2^2^^1^^2  •  (80) 


The  Gaussian  case  will  be  used  here  also  where  05-  =  =  a  and  py  =05-  =  p.  Thus 

X  J  x 2  X  1  X') 


;1  Kx2 


fx,fc,<*l>*2>* 


27T02v/l-p^ 


exp - y- - T  [(Xj-p)2  -  2p(x]-p)(x2-p)  +  (x2-p)2) 

2az(l-p“) 


with  characteristic  function  given  by 


CT2,.  .2 


QxiX2^wl,tJ2^  =  exp| — 2"(w  j  +  2pu\u2  +  w2)  +  JP(C0 1  +  w2)  <8 

The  function  Px  ^  x  (k  j  ,k2 )  is  well-behaved  for  any  pair  of  real  values  (kj  ,k2)-  Its 
Fourier  transform  Qx^^l  ,u>2^  *s  8iven  by 

r  00  r  00 

Qx,x^(a>l’a,2)=  /  I  pXlx^l^2)exP(jwl^l  +  jtJ2^2)d^ld^2 

1  1  J  -00  J  -00  1  a  (8 


=  sinc(cjj/27r)sinc(w2/27r)Qy^-i(cu  j/q,u>2/q)  • 


-4.  • 


Note  that  (4) 


dc71Qx1X2(a,l’<J2) 


=  jp/q 


(Oj,co2=0 


and 


dcojdcop  Qx1x2(cj1’co2) 


_  q2  R*lx2(0) 

tJj,a)2=0  4 


(84) 


(85) 


MULTIPLIER  INPUT  AUTOCORRELATION/CROSS-CORRELATION 

The  multiplier  input  autocorrelation  function  for  non-zero  time  lag  {17T;  17  =£  0}  is 
Rx(r?T).  It  can  be  written  equivalently  as  Rx^x^(0)  which  is  the  cross-correlation  function 

for  two  quantizers  with  equal  quantile  step  size  q  and  zero  time  lag.  The  function  presented 
here  is  valid  only  when  |p| <  1 .  If  p  =  1 ,  the  autocorrelation  function  is  equal  to  the  second 
moment  (76).  Thus, 

Rx1x2(°)=Efxlx21 

-,K 

=  2~2  E[k,k2] 

OO  OO 

=  (12  X  2  klk2Px,x2(kl’k2> 

kj=-°°  k2=-°° 


I 

S|=-°°  s2=— 00 


j1du-lQx1x2^w]’w2) 


(86) 


0)]=27TS] 

CO-)=27TS2 


=  Rx ,  x2(0>  +  ^  Rx ,  x2<®) 


where,  from  (85) 

R„  „  (0)  =  Rtt  ?  (0) 
xlx2  xlx2 


=  po2  +  p2 


(87) 


and.  for  the  Gaussian  case 


oo  oo 


iRxx(0)  =  ^  £  ^  ldii_Jexp(-2(a/q)V(sf-2psiS2+s|))cosj2jr(Sl-s2)/;/q) 


tlX2""  ^  ‘>*2 


oo  oo 


^  ^  till— iexp[-2(£j/q)V(s|+2ps1s2+s^)|cos|2jr(s1-K2>M/q} 


.  q. 

2tt 


Sl=l  S2-1 


+  4 po2  ^  (-l)sexp|-2(a/q)2(7rs)2|  cos^jrsp/qJ 

S=1 


(88) 


+  ^-Mq  ^  exp|-2(a/q)20rs)2  jsin|27TS)i/qJ  . 

s=  \ 


23 


V.  ERROR  STATISTICS 


The  equations  of  the  error  statistics  are  derived  in  this  section.  The  asymptotic  and 
decreasing  terms  are  separated  and  identified.  The  asymptotic  terms  are  found  to  be  inde¬ 
pendent  of  the  p.d.f.  shape.  The  decreasing  terms  are  left  general  in  that  they  contain  the 
forms  Qx  and  Qxx.  This  section  is  intended  to  show  the  derivations  and  results  with  a 

minimum  of  explanation  or  discussion.  Gaussian  forms  of  the  decreasing  terms  are  presented 
in  the  next  section.  Properties  of  the  asymptotic  terms  are  then  discussed  in  the  subsequent 
section. 

If  the  reader  is  concerned  primarily  with  using  the  equations  for  Gaussian  quantizer 
inputs,  he  should  use  only  the  asymptotic  terms  from  this  section  and  the  decreasing  forms 
from  the  next  section. 


PROBABILITIES 

Consider  the  case  of  a  single  multiplier  with  multiplication  coefficient  a,  associated 
parameter  i>,  and  word  lengths  K,  L,  M,  N  fixed  with  K+L  =  M+N.  By  Property  1 ,  the  values 
of  k  which  map  onto  a  particular  value  of  the  error  e  are  given  by  the  equation  k  =  k'+p2^ 
where  the  relation  between  k'  and  e  is  given  by  (58).  In  line  with  the  relaxation  of  restric¬ 
tions  in  the  previous  section,  it  will  be  assumed  k  (and  hence  p)  is  not  bounded.  The  proba¬ 
bility  of  occurrence  of  each  e,  Pe(e)  =  P£(e;  a,u),  can  be  written  in  terms  of  the  probability 

of  occurrence  of  each  k,  Px(k)  which  maps  onto  it.  The  factors  a,  v  and  the  word  lengths 
K,  L.  M.  and  N  all  condition  the  values  computed  of  Pf(e)  for  any  given  p.d.f.  f^(x).  This 
dependence  is  assumed  understood  in  the  following.  Thus. 

OO 

Pe(e)=  ^  Px(k'+p2v)  (89) 

p=— OO 


2  "  ^  Qx(cos)  exp(-jk'cos) 
s=-°° 


=  Pe  +  |Pe(e) 


where 


Pe  =  2~v  ,  (90) 


|Pf(e)  =  2~v  Qx(ojs)  exp(-jk'ws)  ,  (91) 

c=— OO 

s=£0 
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and 


ws  =  2ns/2v  .  (92) 

The  probabilities  can  also  be  derived  for  joint  r.v.’s.  Let  aj  and  a 2  be  the  constant 

coefficients  of  two  separate  multipliers  with  K,  L,  M  and  N  the  same  for  both.  This  results 
in  two  values  of  v,  and  which  are  not  necessarily  equal. 

Consider  the  sets  of  integers  k  j  =  k'j+pj2  1  and  k2  =  k'2+p22  1  where  k'j  ,k'->,»»j  ,n->,  pj 
and  pi  are  integers  with  properties  as  discussed  above.  These  values  of  kj  and  map  onto 
a  pair  of  integers  ej  and  e-i  with  values  determined  by  k'j  and  k2  respectively  by  (58).  The 
probability  of  occurrence  of  each  pair  of  errors  (e  j  ,e->)  is  written  as 
Peje^(ej,ei)  =  Pe  (ej,e2;  aj,ai,i>j,ni).  This  probability  can  be  written  in  terms  of  the 
probability  of  occurrence  Px  x-/kl»k2^  eack  P a‘r  (kpk'O  which  maps  onto  it.  Thus 

OO  OO 


Pe,e-i(el’e2>=  S  2  Px,x^kl+Pl2,;i'k2+P22l,2> 


Pj=_oo  Pt=_oo 


_  2-l/l  "2  ^  £  Ox lX2(«Sl."s2)  times 

c  .  —  ^ c  ~  — OO 


sp*"  s2 


expHk'jw  -jkiu;^) 


=  Peie2  +  |Peie2(el’e2) 


where 


P-  ,  =2  1  2 

e,e2 


^Pe,e^el>e2>=2 


OO  OO 

!  S  S' 

S]  =  -°°  S2  =  _0° 

(sj,s2)it  (0,0) 


0Xlx2K,-ws2)  times 


expHk'jcoS]-jk^coS2) 


o)^  =  2rrs:/2  1  for  i  =  1 ,2  . 

Sj  1 


MEAN  ERROR 

Through  (89)  the  mean  error  n€  can  be  expressed  as 
Me  =  E[e] 

=  2-M-t>  gjej  (97) 


=  2~M-r>  ^  e  Pe(e) 
e 

=  +  lfie 


where 


e 


2  M  v-\  for  TCR  (Two’s-Complement  Roundoff) 
2-M-l(2~r'_i)  for  jcc  (Two’s-Complement  Chopping) 
and 


(98) 


iu€  =  2'M-"  ^  e  (Pc(e)  . 
e 


(99) 


The  sum  (98)  is  evaluated  through  use  of  (62)  and  (63).  For  convenience  the  limits  on  e  are 
not  shown  because  they  depend  on  whether  roundoff  or  chopping  is  used. 

The  mapping  from  k'  to  e  is  one-to-one.  Hence,  the  error  can  be  written  as  a  function 
of  k';  namely  e  =  e(k').  Also,  a  summation  over  all  possible  values  of  e  is  equivalent  to  a  sum¬ 
mation  over  all  possible  values  of  k'.  Thus,  through  (9 1 )  the  term  can  be  written  as 


ine  =  2'U~v  J  e(k')  ^  2~v  Qx(ws)  exp(-jk'ws) 

e  s=-0° 

s=£0 


OO 


_  2~M-2t' 


I 


«--00 

s#0 


QxK>Ds 


(100) 

(contd) 
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where 


=  2-M-2H-l 


£  Re{QxK)Ds} 

S=1 


2"-l 

Ds=  ^  e(k')  exp(-jk'cos)  . 
k'=0 

The  quantity  Ds  is  evaluated  in  the  Appendix  and  is  given  as 
h"-1  for  s  =  0  mod  2V 


Ds  = 


for  TCR,  and 


Ds=  < 


.  exp((sgn  a)jco  X)-l 
(— 1  )s  2v~l - ~~~r~  ~ — ; -  otherwise 


cos(cosX)-l 


2^1(1  -2*')  for  s  =  0  mod  2V 


2^-1 


exp((sgn  a)jcJsX)-l 
cos^A)-! 


otherwise 


(101) 


(102) 


(103) 


for  TCC.  Note  that  Ds  can  be  a  complex  quantity  and  that  it  is  a  function  of  the  parameter  X 
which  satisfies  the  relation 


1  =  |C'|X  mod  2V 

where  8'  is  in  turn  related  to  the  value  of  the  coefficient  a  through  (47)— (52). 


(104) 


Re(0) 

The  second  moment  of  the  multiplication  error  for  one  multiplier  is  denoted  by 
Re(0).  Thus 

Re(0)  =  E[e2] 

=  2~2M— 2l/  E[e2] 


=  2-2M-2k  ^  e2Pe(e)  (105) 

e  (contd) 
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=  Rg(0)  + |Rg(0) 


where 


R€(0)  =  2~2M'3,;  ^ 


-7-2M 


[2-3-2~v +  2~^p] 


and 


for  TCR 

for  TCC 


*Rg(0)  =  2~2M'2t'  J  e2|Pe(e) 
e 


=  2~2M-2 v  e2(k')  ^  2  v  Qx(ios)  expt  -jk'a\.) 
e  s=~°° 

S9 to 


,2-2M-3>m-I  £  Re{<3x(Us)Fs). 

S=1 


The  quantity  Fs  defined  by 


Fs=  S  e2<k')  exp(-jk'cos) 
e 


is  evaluated  in  the  Appendix  and  is  given  by 


±(23p  +  2’21’)  for 

s  =  0  mod  2‘ 

■jV 

f-l)s  2 

^  '  l-cos(u>  X) 

otherwise. 

for  TCR,  and 


(106) 


(107) 


(108) 


(109) 
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i 


F 


s 


±(2-23v -3-22v +  2V)  for  s  =  0  mod  2V 


2~t'~*  [exp((sgn  a)ju?sX)-l  ]  +2P 
l-cos(cos\) 


otherwise 


(110) 


for  TCC.  It  can  be  complex  and  is  dependent  on  the  coefficient  value  through  the  parame¬ 
ter  X  defined  by  (104)  which  is  the  same  relation  as  used  for  Dg. 


Rxe  for  p  =  1 .0 

The  joint  moment  of  a  r.v.  x  and  the  corresponding  error  e  for  the  same  multiplier 
is  denoted  by  Rxe  where 

Rxe  =  E[xe) 

=  q2-M-i/E[ke]  (!!!) 


and 


OO 

E[kel  =  ^  ke(k)Px(k) 
k=-°° 


OO  2v-\ 

■1  I  (k'+P2p)  e(k'+p2t')  Px(k'+P2") 

p=-°°  k'=0 

2v-\  00 

=  J  e(k,)  S  (k'+p2^)  Px(k'+p2,;) 
k'=0  p=-°° 


2P-1 

=  ^  e(k')  ^  2-,/(-j)|^Qx(cj)|  exp(-jk'tos) 
k'=0  s=  ~ 


s=- 


-2-%D0/,-j2- J  j£Qx(U)}^  Ds 


s=- 
s=£0 


(112) 
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The  third  step  in  (1 12)  is  possible  since  e(k')  =  e(k'+p2f/)  and  since  the  mapping  from  k'  to  e 
is  one-to-one.  The  last  step  is  obtained  through  use  of  (69).  Thus, 


^xe  ^xe  +  ^xe 


(113) 


where 


Rxe-P2'M-2>’D0 

forXCR 

m2-M-i^1(i_2»')  for  TCC, 
and 

OO 

£  {£<3xto>}  Ds. 

s=-°o  S 

s#0 


(114) 


(115) 


Note  that  this  is  the  same  joint  moment  as  for  the  r.v.  Xj  associated  with  multiplier  a  j  and 
the  error  e2  associated  with  multiplier  a2  when  p  =  1 .0.  In  this  case  it  is  necessary  to  make 
the  assignments  v  =  i>2  and  s  =  s2  so  that  ws  =  2rrs2/2  -. 


Rxle2  f°r  <  10 

Let  x  j  be  the  input  for  the  multiplier  with  coefficient  aj  and  e2  be  the  error  for  the 
multiplier  with  coefficient  a2-  Furthermore,  let  |p|  <  1.  The  joint  moment  of  the  r.v.’s  xj 
and  e2  is  denoted  by  Rx  e  where 


oo 


oo 


E[k,e2]  = 


l 

; 


2  i  kie2(k2)px  x2(ki.k2> 

kj=-°°  k2=-°° 


oo  oo  2^-1 


III  kl  e2(k2+P22  2)Px,x2(kl>k2+P22  2) 

kj=-°°  p2=-°°  k2=0 


2VU 


oo  oo 


S  e2<k2>  2  2  kl  PXjx2(kl>k2+P22  2) 


(117) 


k2=0  kj=-°°  p2= 


2  2-l 

I 

k2=0 


e2(k2)2  ^  times 


oo  oo 


I  2« 

Sj  =  -oo  s  =_<*, 


dco,  QXjX2(wl>w2) 


COj=27TSj 

"2^8, 


exp(-jk2coS2> 


oo  OO 


-i  "2  2  I 

Sj=-°°  S2=-< 


dco,  Qx1x2(ail’co2) 


D, 


co  j=2jtsj 

co2=cos. 


s2 


-}2 


n~v2 


OO  OO 

2  2 

Sj=-°°  s2  =  -°° 
(S],S2)  ^(0,0) 


doo^XjX^l’^) 


Dc 


COj=27TSj 

w2=<js2 
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i*', 


where 


v-, 

2  2-l 


Thus 


where 


and 


Ds  =  Z  e2(k'2)exp(-jk2ujSi)  . 
k'2=0 

(118) 

Rx,e2  =  Rxj£2  +  iRXje2 

(119) 

Rxje2-#i2  ~[Ds2}s2  =  0 

(120) 

p2  2  forTCR 

-M-v-i-l  v~t 

p2  2  (1-2  2)  for  TCC, 

(121) 

;Rxle2  = 

OO  OO 

S  ]  =  — 00  Si  =  -°° 

(s,,s2)#(0,0) 

£:QX|X,'"I"2>  Ds,_ 

1  cj]=27rsj  - 

w2=a;s: 

(122) 

The  joint  moment  of  the  r.v.  x  and  the  error  e  for  the  same  multiplier  when  |p|  <  1 .0  is 
obtained  by  simply  letting  aj  =  a 2  in  the  above. 

Re,e2f°r  P=  10 

Let  ej  and  e2  be  the  multiplication  errors  for  the  multipliers  with  coefficients  aj  and 

a2,  respectively.  Also,  let  p  =  1.0  and  >v~,.  The  joint  moment  of  the  r.v.’s  ej  and  e2  is 

denoted  by  R,  ,  where 
ele2 

Reie2  =  E[eie2] 


=  2 


-2M-»'j-i'2 


E[e j e2 ]  . 


(123) 
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The  number  of  distinct  states  of  the  pair  (e|,e2)  is  determined,  in  this  case,  by  »»j,  since 
>  vi.  The  probability  of  occurrence  of  each  state  is  the  probability  of  occurrence  of  ej 
alone.  Also,  since  p  =  1.0,  the  error  value  e^  is  uniquely  determined  by  the  error  value  e  j . 
Thus, 


V^'2M  ^  J  e,  e2(k'1(e1))Pei(e,) 
el 

~  +  j e 2 


(124) 


where,  from  (89) 


Re\€2 


-2M-2c  i  — v-> 


1  V2  £  el  e2(k'|(e j )) 


2^-1 

=  2-“N1-“^  1  _l"2  e,(k'| ) e2(k  j ) 

kj=0 


and 


IK 


le2 


-2M-2i'|-i't+ I  \' 

=  2  1  ~  /  el  k j (e j  ))  times 


J  Re  {Qx(u;s  i  >  exP< — J k  !  * e  I  K, 


sl  =  l 


Vi 

2  ‘-1 


-2M-2(',-r".  +  l  V  , 

=  2  -  2L,  e  j  ( k  j  )e -» ( k  j )  times 


k',=0 


(125) 


( I  2(>) 


^  Re  Qx(coSj )  exp(-jk'|UJS)  )J 
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Example  plots  of  ej  versus  e^  appear  in  [61 .  Unfortunately,  attempts  at  arriving  at  more 

suitable  closed  forms  for  it-  .  and  |R_  -  have  been  unsuccessful  except  for  the  special 

clc2  ele2 

case  ej  =  ei  already  derived.  (For  =  e^,  see  the  equation  for  Rc(0).)  This  is  the  reason 
for  reverting  to  sums  which  are  more  easily  computed  using  the  index  k'|  instead  of  ej . 


Rfje2  for  Ipl  <  10 

Let  e  j  and  e->  be  the  multiplication  errors  for  the  multipliers  with  coefficient  values 
a|  and  a->.  respectively.  Also,  let  jp|  <  1  as  for  RX|6V  The  joint  moment  of  the  r.v.’s  ej 

and  e->  is  Rf  where 

Re,e2  =  Etele2l 


E[e,e:l 


-2M-v,-i'1  V" 

=  -  ■  Z  Z  ele2  Pe,e2  <el  e2> 


e  i  ei 


R-  -  +  iR,  , 

€  |  €  -l  T  €  1  €  i 


where,  from  (93) 


-*  -2M-2i>|  -2ri  V  V 

V:'-  ’  -  2  1 


(127) 


e,  e-, 


for  TCR 


2  1  -  (1-2  1  >(  1-2  ->  for  TCC 


and 


(128) 


Through  (95)  the  term  J.R 


can  be  written  as 


■\e2 


M-y  i  XT 

|Reie2=2  -  Z  Z  el(k'l)e2(k2)  times 

el  e2 


OO  CO 


S  S  ^  t/2Qx1x2^s1-ws2)exPHk'1roSi-jk'2wS2) 

Sj  =  -°o  St=-°° 

(s,,s2)^(0,0)  (130) 


OO  OO 


=  2-2(MWi+»2,  £  £  QX]Jl7(toS|,Us2IDS|DSj 


Sj  =  -oo  St  =  -°o 

(Sj,S2)#(0,0) 


The  autocorrelation  of  the  error  for  one  multiplier  can  be  obtained  from  Rf(0)  for  p  =  1 .0  and 
(130)  for  |p|  <  1.0  by  letting  aj  =  a->. 


VARIANCES  AND  COVARIANCES 

The  variance  of  the  multiplication  error  Ce(0)  is  defined  as 

Ce(0)  =  Re(0)-p2 

=  Ce(0)  +  lCe(0)  (131) 

where,  from  (97)  and  (105), 

Ce(0)=  Re(0)-  M2  (132) 


for  TCR  and  TCC,  and 

lCe(0)  =  iRe(0)-2peipe-ipJ  .  (133) 
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The  function  C'x  £i  is  the  covariance  of  the  multiplier  input  Xj  (for  the  multiplier 
with  coefficient  a|  1  with  the  error  t-»  (for  the  multiplier  with  coefficient  a->).  Similarly,  the 
function  C'Xf  is  the  covariance  of  a  multiplier  input  with  the  corresponding  error.  Both  can 
be  defined  as  the  following  (dropping  the  subscripts  in  the  first  case  for  brevity): 

^  xe  -  Rxe  ~^\^e 


=  txe  +  iCxe 


where,  from  (72).  ( 1 13)  and  (1  19), 


^  xi  ~  R\e  ^x^t 


for  t  C  R  and  K  (  .  and 


iC'X(:  =  f  Rxt  ~  ~  U3f 

rite  covariance  C'c  t  ^  of  the  multiplication  error  e  j  (for  the  multiplier  with  coeffi¬ 
cient  a | )  with  the  multiplication  error  e-,  (for  the  multiplier  with  coefficient  a->)  and 
Ipl  <  1 .0  is  defined  as 

C",  c  =  R,  ,  -  u,  (1 3' 

6  1 1  -»  6  |  6  -»  't  j  i 


=  c*  +  , 

e  j  e ->  e  |  e-i 


where,  from  (97 )  and  (127) 


C,  .  =  R,  ,  -  juc 
e |6-i  e  |  e  -i  f|  fi 


for  TCR  and  TCC,  and 


iCe,e2  = 


1  fc  l  fc2 


When  p  -  1.0.  the  asymptotic  portion  Ce  £  is  generally  not  equal  to  zero.  However,  the 

i  2 

decreasing  portion  fC,  ,  has  the  same  form  as  ( 1 39). 

j  ^ 


Tlie  following  comments  regard  the  asymptotic  correlation  coefficient  (ACT)  pc 
defined  for  p  =  1 .0  as 

'  "401 

(A  similar  form  of  the  ACT  was  examined  by  Girard  K>|  and  Parker  and  Girard  17].  They 
assumed  a  zero-mean  error  which  is  not  the  case  for  TCR.)  f  irst  consider  the  case  when 
P|  =  p->  =  p.  f  or  a  given  v  the  ACT  values  vary  from  a  maximum  of  1 .0  to  values  which  can 

become  quite  small  but  which  are  non-zero.  Table  1  shows  the  computed  ACT  values  which 
result  when  v  =  5  and  for  two  coefficient  values  Gj .  Note  that  the  same  ACT'  values  result 

for  both  cases  but  they  are  permuted  with  respect  to  the  values  of  C'-t.  This  permutation 
property  appears  to  be  a  general  one  for  all  odd  values  of  C'| .  Also  given  in  the  table  are  the 

computed  mean  and  standard  deviation  of  these  ACT  values  which  are  therefore  the  same 
for  all  odd  values  k'[ .  Table  2  shows  the  computed  mean  and  standard  deviation  of  the  ACC 

values  as  a  function  of  p  for  p  up  to  10.  The  decrease  in  the  standard  deviation  value  as  v 
increases  results  front  the  introduction  of  more  ACC  values  which  are  closer  to  zero. 

The  next  case  to  consider  is  when  p  |  >  p-i.  Table  3  shows  some  example  ACC  values 

which  illustrate  the  similarities  and  differences  resulting  from  the  different  coefficient  values. 
Ihe  common  behavior  is  for  the  ACT'  values  to  start  at  some  initial  value  (for  Pi  =  5  in  this 

case),  then  eventually  decrease  by  a  factor  approaching  2.0  for  each  integer  increase  in  Pj . 

The  differences  in  the  ACC  values  can  be  seen  in  the  four  cases  which  are  shown.  The  cases 
were  chosen  according  to  the  signs  of  the  ACT  values  for  Pj  =  5  and  P]  =6.  respectively. 

f  our  possibilities  (++.  +-  .  +  and  -  - )  are  represented  here.  For  pj  >  b,  no  polarity  changes 

are  indicated  and  the  decrease  in  absolute  value  of  the  ACT  by  factors  close  to  2.0  takes 
over. 
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VI.  GAUSSIAN  FORMS 


In  this  section.  Gaussian  forms  of  the  decreasing  terms  of  the  error  statistics  are 
written  out.  They  were  obtained  by  substituting  the  Gaussian  forms  of  Qx  and  Qxx  (from 

(08)  and  (83))  into  the  decreasing  terms  defined  in  the  last  section.  In  the  interest  ot  a 
simplified  presentation,  a  shorthand  notation  is  employed  witli  the  following  definitions: 

A  =  -(o/q)~/ 2 


B  =nl  q 


G 


( -  T)s  for  TCR 
1  for  TCC 


H 


1  for  TCR 
1  -2V  for  TCC 


(141  ) 


V  =  l/(cos(cesX)  -  1 ) 
W  =  sine  (  ojs/27t  ) 
cos=  2ns/21'  . 


Subscripts  are  used  with  G,  H.  V.  W  and  cos  whenever  more  than  one  coefficient  is  involved 
in  the  formula  of  the  statistic.  These  have  the  form 


C: 


H; 


(-1  )S‘  for  TCR 
1  for  TCC 
1  for  TCR 
\-2V{  for  TCC 


Vj  =  l/(cos(cus.Xj)- 1 ) 


W:  =  jinctcu  /27r) 

1  Sj 


to 


v- 

2TTSj/2 


(142) 


for  i  =  1.2.  There  is  a  special  term  coi  used  for  J.R 


X|Ci 


for  |p|  <  1 .0. 


In  t his  case 


It  is  assumed  that  the  coefficients  have  fixed  values.  Thus,  whenever  the  parameter  k' 
appears,  a  specific  transformation  from  k'  to  e  is  implied  according  to  (58).  The  parameter  X 
is  the  coefficient  related  value  obtained  from  (104). 


PROBABILITIES 


iPf(e)=2*  v  ^  W exp(A<o"> cos(a>sl  B -k’) ) 
s=l 


(1431 


*Pf,e2,el-e2»=  2 


I 

s,  =  l 


W,  exp(Acj-  )  cos(cj„ 


I B - k'j  | )  +  2,  W2  exp(Aw^)  Cos(u)SlIB-k’;|  l 
Sl=l 


OO  CO 


X  Iwiw:T' 

S  |  =  1  S->=1  ^ 


+  exp 


exp  |  A(cjJ  +  2p  cjs  )  cos(cjs  ^  I  B-k'j  1  +  cjs>|  B  -  kM  ) 


(144) 


{A(‘js|-:R<js,<js:  +  wsv}  eos(«S||B-k’, 

J 

MEAN  ERROR 


TCR  and  TCC 

oo 

*P£  =  2-M“'  £  G  VW  explAto")  [cos(ws|B  +  k  sgna)  )  -cos(u>sB>] 
s=l 

s  ^  0  mod  2" 


(145i 


Re(0> 

TCR 

oo 

IR,(0)  = -2':M':v+I  y  (-l)s  VWexpfAurHosft^B)  .  (146) 

c  * — «  S  3 

s=I 

s£0  mod  2" 

TCC 

oo 

iRf(0)  =  -2~2M-2v+I  ^  VW  expfAcjj  ) |2,,_1  cos(cjs| B  +  X  sgnal )  +  ( I  -2P~' )  cos(ojsB)1  . 

s=l  L  J 

s^O  mod  2r  (147) 
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/ 


Rxe  for  p  =  1 .0 


TCR  and  ICC 

W>  oo 

i-K-M-i/ 1|  J-exp(Aoj-)LOS(cjs/2)sin(cjsB)  +  ^  exp(Au>“)  times  (148) 


‘Kxt  = 


W 


s=l 
s=0  mod  2U 


S=I 

s  ^0  mod  Z” 


+  —  cos(o)  J2) 


VG  sin(io..|  B  +  X  sgna]  )  -  sin(tosB) 


■  GBVwj  cos(a>s|  B  +  X  sgna) )  -coslu^B) 


Rx  t  for  |p|  <  1.0 


TCR  and  TCC 

oo 

r,  -  V  if  I  I 

iRX|e-,_-  ‘  Z.  G2W2  exP(Aa,s-,)  |^-AV2PCl>s-,  |sin,‘Js-.lB  +  ^2SBna2l,-s'n(tos-,B)} 


Sl=l 
“  Pt 

s->^0  mod  2  ~ 


+  BVi  cos(uSiIB  +  X-,  sgna-i| )  -  costw^Blj 


-K-M-e-,,,  V  ■>  (_|,sl 

+  -  '  Hi  £  exp(Aco]  )  - — - —  sin(co|  B) 

>< 


OO  OO 


-K-M-Vi  V  V  s, 

~  2*  /,1'D  Gi  exp  A(co][  +  2pu>i  cas^  +  w"^)  times 


(149) 


s  |  =  I  s->=  I 

v-i 

St  =f0  mod  2  * 


WtV, 


Jsin(  B(<j|  +  ojSi)  +  uSiXi  sgna-> )  -  sin(B|u>  |  +  I  )J 


OO  OO 


+  2 


-K-M-v-, 


^  2  (-1 1  *  Gt  exp  j  A(u>  |  -  2pui  |  u)Sl  +  )  j  times 

sl  =  !  s2=l 

*>-> 

sj^O  mod  2  - 


W-.V-, 


—  jsin(B|<j|  -  ojSi)  -cj^Xt  sgiian)  -sin(B[oj|  -  coSnl  )J 
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Reie2forP  =  10 


TCR  and  TCC  and  v^>V2 


IR, 


el£2 


v, 

2  '-1 

I 

kj=0 


e  j  ( ) 


i 

Lsr' 


Wj  exp(Au>^)cos(toSi  (B  -  k'j  I ) 


(150) 


Rele2f°r  'p*<  10 
TCR  and  TCC 


IR 


ele2 


+  2 


-2M-p  .  1  V"  t  r  1 

2  1  *■  H2  ^  GjV|Wj  exp(Aw^  )  jcos(a)SjlB  +  Xj  sgnaj  1 )  -  cos(cjs^B)J 

si=1 

Sj  P  0  mod  2  * 

OO 

-OM-Vi-i/v-l  V  ■>  r 

Hj  2  g2V2W2  exptAco^)  coskos,lB  - 

s2=l  *  2  L 

s->^0  mod  2  * 


(  +  X-i  sgna-,  | )  -  cos! u> 


s:B>] 


+  2 


ST  f  -t  \ 

2  2GlG-ViV2WiW,expjA(u,sr  +  2pcos  cos^  +  u>“^)  Mimes 

s,=l  s2=l  '  *  1  *  2 

S|  ^0  mod  21*1  - 
v-i 

s2#  0  mod  2  -  (I5D 

^cos(B(ojSi  +  coSll  +  <0^  X)  s8nal  +  ws.,X2  sgnaT )  -  cost  B[u>s^  +  (0^1  +  wS)X|  sgna]) 

-  cos(B(  wSj  +  10^)  +  X2  sgna2)  +  cos(B|u>S)  +  co^l  )J 


•2M-i'|-i>i-l  ^  V  |  ■>  -1  ) 

2  2  gIg2v1  V2  wl  w2  exp  (A(u,si  '’pu’s1tJs-.  +wsi  M 

st=l  s2=l  *  ‘ 

S|  ^0  mod  2  ' 
st^O  mod  2V“ 


times 


^cos(B(wS)  -oo^l  +  wS)  X  j  sgnat  -c;s,X2  sgna-,) -cost  B|<oSi  -ws-)]  +  co^  X(  sgnat ) 


-cos(BlcoS!  -w^I  -toS2X2 


i2)  +  cos(B|coS) 


41 


VII.  DISCUSSION 


It  is  appropriate  at  this  point  to  compare  the  results  derived  in  this  report  with  the 
error  models  presently  used  in  digital  filter  design.  Obviously,  simplification  of  the  derived 
equations  are  not  possible  when  a/q  is  small  enough  so  that  the  decreasing  form  of  each 
statistic  cannot  be  ignored.  For  the  purpose  of  this  discussion,  a/q  will  be  assumed  large 
enough  so  that  the  decreasing  forms  are  negligible. 

Assume  a  filter  structure  where  the  coefficient  values  are  fixed.  Associated  with  each 
coefficient  value  is  the  parameter  v  which  is  the  effective  word  length  of  the  associated 
multiplication  error.  The  errors  can  take  on  each  of  2 v  discrete  values  which  are  uniformly 

distributed.  The  probability  of  occurrence  of  each  value  is  (almost)  equal  to  2~v .  The  mean 
— ► 

error,  ne,  is 

— y 

»e  = 


1- M-n-l  for  XCR 

2- M-1(2-j'_i)  forTCC 


(152) 


The  error  variance,  Ce(0),  is 
-  i-2M 

Ce(0)  =  ^T2~(1-2  } 

An  idealized  continuous  uniformly  distributed  error  ec  is  often  used  in  filter  design.  This 
error  has  the  following  properties: 


(153) 


Pe 

Cc 


0  for  TCR 
-2~M_1  forTCC  , 


(154) 


and 


i-2  M 

c<»=iF  <IS5> 

Note  that,  in  comparing  ( 1  54)  with  (152)  and  (155)  with  (153),  the  idealized  continuous 
model  can  only  be  assumed  if  v  is  large  enough.  Thus,  for  example,  if  the  word  lengths  K.  L, 
M  and  N  are  equal,  the  coefficient  values  ±1/2  (or  v  =  1 )  would  yield  the  most  difference 
between  the  discrete  and  continuous  model.  The  difference  becomes  less  for  coefficient 
values  of  ±1/4  and  ±3/4  (or  v  =  2),  and  so  on  for  larger  values  of  v. 

Other  results  are  the  following.  The  covariance  Cxe  is  zero  in  agreement  with  the 

continuous  model.  This  means  the  multiplier  input  is  uncorrelated  with  the  error.  The  error 
covariances  are : 


C£(t?T)  = 


C£(0)  for  r?  =  0 
0  for  r?  #  0 


(15b) 


i 

i 

i 

i 

i 


and 


_  CeiC2  forp-LO 

jCt  ~  * ' ^7) 

l  0  for  |p|  <  1.0  . 

These  latter  results  imply  that  these  contributions  to  the  FIR  filter  variance  and  output 
spectrum  are  white  noise  in  nature. 

The  consequences  of  small  values  of  o/q  on  the  statistics  will  be  explored  in  a  sub¬ 
sequent  report.  Particular  attention  will  be  paid  to  those  ranges  of  values  of  o/q  over  which 
the  asymptotic  model  can  be  assumed . 


Mean  of  p 


le2 


Standard  Deviation 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


1.0 

0.6 

0.33333 

0.17647 

0.09091 

0.04615 

0.02326 

0.01167 

0.00585 

0.00293 


0.0 

0.4 

0.4762 

0.4216 

0.3566 

0.2737 

0.2052 

0.1497 

0.1086 

0.0777 


Table  2.  Computed  mean  and  standard  deviation  of  P€^e^  for 
TCR  as  a  function  of  j>j  =  v2  =  v. 


v\ 

8  j  =  17 
(++) 

1  II 

Cj  =  19 
(-+) 

8’,  =  15 
(-) 

5 

0.249267 

0.354839 

-0.020528 

-0.067449 

6 

0.124588 

-0.104067 

0.083547 

-0.127519 

7 

0.062288 

-0.052029 

0.041770 

-0.063754 

8 

0.031143 

-0.026014 

0.020884 

-0.031876 

9 

0.015572 

-0.013007 

0.010442 

-0.015938 

10 

0.007786 

-0.006503 

0.005221 

-0.007969 

NOTE:  j>2  =  5,8'2  =  1 


Table  3.  Examples  of  the  behavior  of  p£  e  f°r  TCR  when  P| 
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Block  diagram  representation  of  an  FIR  filter  which  incorporates 
quantization  errors. 


SIGN 

BIT 


MULTIPLIER  INPUT 
x 

A _ 


K  BITS- 


"V 

4- 


COEFFICIENT 

a 

_ A _ 


A 


■L  BITS  - 


BIT  POSITIONS 
EQUAL  TO  ZERO 


V  BITS 


00 


8  BITS 


00  • • •  0 


■M  8ITS- 


-N  BITS- 


- v - 

y 

PORTION  LEFT  OVER  AFTER 
ROUNDOFF  OR  CHOPPING 


-A- 


AA 


PORTION  DISCARDED  IN  THE 
ROUNDOFF  OR  CHOPPING 
OPERATION 


Figure  2.  Relationships  among  the  various  word  lengths. 
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APPENDIX  A 


Evaluation  of  Ds  proceeds  as  follows.  Consider  the  TCC  case  first  where  a  0.  Then 
by  using  (58)  in  ( 101 ). 

2*'—  1 

Ds  =  ^  e+(k')  exp(-jk'ccs) 
k'=0 

2"-l 

=  -  ^  [k'C'l  it>  exp(-jk'cjs)  (A-l) 

k'=0 

for  tus=27rs/2t'  where  s  =  0,±  1  ,±2.  .  .  . ,  and  C'  is  the  coefficient  dependent  value  which  is 
derived  from  the  value  of  a  through  (47),  (51).  Let  the  variable  j3  be  defined  by  relation 

(3  =  k'C'  mod  2V  (A-2) 

which,  since  C'  is  an  odd  positive  integer,  implies  the  relation 

k'=0Amod21'  (A-3) 

where  the  odd  constant  integer  AG[ 0,2")  is  the  unique  inverse  of  C'  which  satisfies  the  relation 

1  snmod  2"  .  (A-4) 

Since  the  relationship  between  k'  and  P  is  one-to-one, 

2v-\ 

Ds=-^  P  exp(-joos[0A]  2^) 

13=0 

j 

=  -/  P  exp(-jcos|3X)  .  (A-5) 

P=0 

Whenever  s  =  0  mod  2V  this  form  simplifies  to 
2"- 1 

Ds«-  £  0 
0=0 


=  2V-I( 1-2^)  . 


(A-6) 


This  form  can  be  evaluated  for  other  values  of  s  by  writing  it  as 


2u-\ 


D  =  y  -S--5S-1 
s  zl,  jco^.  dX  1 
0=0 


exp(-jcos 


1  d 
jcos  dX 


2v-\ 

I 

0=0 


exp(-ju;,.0X) 


( A-7 ) 


The  summation  is  over  a  power  series  in  exp(-jtosX)  which  is  easily  evaluated.  By  then 
taking  the  derivative,  the  result  is 


D 


s 


^1  exp(jcusX)-l 
eos(cosX)-l 


( A-8) 


for  s  p  0  mod  2". 

The  computation  of  Ds  for  a  <  0  makes  use  of  the  relationship  of  e+(k)  to  e_(k)  as 
called  out  in  Property  2.  Thus,  tor  a_  where  a_  =  -a+. 


I3s(a_) 


^  e_(k')  exp <  — jk'cos ) 
k'=0 


=  0  •  exp(O)  +  ^  e+(2^-k')  exp(-jk'cos) 
k'  =  l 

e+(k')  exp(jk'cos) 

k'=0 

=  D*(a+)  (A-9) 


where  *  denotes  complex  conjugate.  As  a  result,  the  complete  expression  of  Ds  for  the  TCC 
case,  which  takes  into  account  the  sign  of  the  coefficient  value,  is  written  as 


Ds=  < 


2*^~ 1  <  1 — 2*')  for  s  =  0mod  2" 

j  exp((sgn  a)jccsX)-l 


cos(cosX)-l 


otherwise 


(A-10) 


i 


The  parameter  X  is  then  determined  from  ( A-4)  whereby  1 1' |  is  used  in  those  eases  where 
a  <  0. 

This  complex  conjugate  relationship  holds  for  Ds  in  the  TCR  case  and  1  s  in  the  TCR 

and  TCC  cases.  Hence,  although  the  bulk  of  the  remaining  derivations  are  carried  out  for 
a  >  0,  the  final  results  in  each  case  will  be  stated  so  as  to  include  the  effect  of  the  sign  of 
the  coefficient  value. 

Evaluation  of  LT  for  the  TCR  case  follows  the  same  procedure  as  for  the  TCC  case. 
Substitution  of  (58)  into  ( lOi )  for  a  >0  results  in  the  following  form. 

Ds=iDs  +  :Ds  (A-ll) 


where 


:"-i 

lDs  =-  X  IkY  +  ^-'l^expt-jk'cOs) 

k'=0 

and 

2"-l 

sDs  =  2i;~1  ^  exp(-jk'cos) 

k'=0 

f  s-"-'  for  s  =  0  mod  2V 


^  0  otherwise  . 

For  the  evaluation  of  ]  Ds  let  the  variable  (3  be  defined  as 


( A- 12) 


t  A- 1 3 ) 


(3  =  (k'E'  +  2V~ 1 )  mod  2V  ( A- 1 4 ) 

which,  since  C'  is  a  constant  odd  integer,  implies  the  relation 

k'  =  (0X  +  2('~* )  mod  2^  ( A- 1  5 ) 

where  the  odd  constant  integer  X  is  the  unique  inverse  of  C'  which  satisfies  the  same  relation 
(A-4)  as  for  the  TCC  case.  With  this  change  in  notation 

2^-1 

1  Ds  =  —  ^  (3  exp  /-jcos[/3X  +  2"-*  i  (A- lb) 

,3=0  '  “  '  (contd) 
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I 


=  -  ^  0  exp  (— jcos[/3X  +  2V~ ' ) ) 

p=o 

~<V_  j 

=  -(-l)s  ^  0expHcos0X) 

0=0 

which  is  exactly  the  same  form  as  Ds  in  the  TCC  case  except  tor  the  tactor  (-1  )s.  Thus,  tor 
TCR,  and  taking  into  account  the  sign  of  the  coefficient  value, 


2V  '  for  s  =  0  mod  2V 

,  „  ,  exp((sgn  a)jcu  X)-l 

(-l)b2y  1  - - - rr — ; -  ,  otherwise 

cos(cosX)-I 


( A- 1  7 ) 


for s  =  0,±1 ,±2,  ...  . 

Evaluation  of  Fs  proceeds  as  follows.  Consider  the  TCC  case  first  where  a  >  0.  Then 
by  using  (58)  in  (108) 


F$  =  ^  e-(k')  exp(-jk'cos) 
k'=0 


( A- 1 8 ) 


=  ^  |  -[k'C'l  w  -  exp(-jk'cus) 

k'=0 

for  co$  =  2tts/2v  where  s  =  0,±  1  .±2 . and  6’  is  derived  from  the  value  of  a  through  (47), 

(51).  Let  the  variables  0  and  X  be  defined  by  the  relations  (A-2)-(A-4)  as  for  Ds.  Then 


Fs=  ^  02exp(-jws[0X]  J 
(3=0 


=  ^  0~  exp(-jcus0X)  . 


(A-19) 


Whenever  s  =  0  mod  2V  this  form  simplifies  to 


0=0 


=  -L  /  s  .  j  2r'  +  -ii>  > 

6  V"  -  "  '  • 

This  form  can  be  evaluated  for  other  values  of  s  by  writing  it  as 


(A-20) 


Fs  =  ~ 


_1^  d^ 
co-  dA“ 


it* 


^  exp(-ju>s0A) 
0=0 


(A-21) 


and  evaluating  the  second  derivative  of  the  power  series  in  exp(-jcosA).  The  result  is  written 
for  TCC  as 


2>~1  [expUcjsX)-l  1+2*1 
l-cos(cosA) 


(A-22) 


for  s  ^  0  mod  2V .  The  complete  expression  which  takes  into  account  the  sign  of  the  coeffi¬ 
cient  value  is  written  as 


^-(2*23i;-  3*2-"  +  2V)  for  s  =  0  mod  2V 


Fs  =  < 


2~v  1  [exp((sgn  aljco^X)-!  ] +2^ 


l-cos(cosA) 


,  otherwise 


(A-23) 


for  s  =  0.±  1  .±2 . 

Substitution  of  (58)  into  ( 108)  for  a  >  0  results  in  the  following  form  for  TCR- 
Fs=  lFs  +  2Fs  +  3Fs  (A-24) 

where 

i 

1  Fs  =  2  |  -  [ k'c'  +  2U~]  1  ^1“  exp(-jk'cos)  ,  (A-25) 

k'=0 
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Am 


I 


where 


2Fs  =  -  2  ^(k'C'  +  Z^'l^expHk'u;,) 
k'=0 


2"- 1 

3FS  =  ^  (2*'-1 )”  exp(-jk'cjs) 

k'=0 


2^"  -  for  s  =  0  mod  lv 


,  otherwise  . 


Let  the  variable  0  be  defined  by  the  relation  (A- 14)  as  for  Ds.  Then 


(A-26) 


(A-27) 


— _ i 

1  Fs  =  J  ^  eXp  ("K  +  2V~]  1  ) 

0=0  “  ' 

2*'-l 

=  (-l)s  ^  02  exp(-ju;s0X) 

0=0 


( A-28) 


which  is  exactly  the  same  form  as  Fs  in  the  TCC  case  except  for  the  factor  (-1  )s.  Thus 
~(2-23v -3-221' +  2")  for  s  =  0  mod  2V 


22*'-1  [expfjco^X)-  1  ]  +2V 


( A-21)) 


1  -cos(ccsX) 


.otherwise  . 


The  factor  ->F  can  be  written  from  (A- 1  2)  as 


2fs  =  2**  ,DS 


(1-21')  for  s  =  0  mod  2V 


(-ir2 


s22^i  exp(jg)sX)-  1 

cos(cusX)-l 


,  otherwise  . 


(A-30) 


The  final  result  is 


(23l'  +  2-2*'J 

for  s  =  0  mod  2V 

■>v 

(-1)S  - - = - -  , 

otherwise 

1  -cos(co  X) 

(A-31 ) 


tor  s  =  0,±1  ,  ±2 . Note  that  this  form  is  real  and.  hence,  there  is  no  need  to  take  into 

account  the  sign  of  the  coefficient  value. 
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GLOSSARY 


ACRONYMS 

ACC 

Asymptotic  correlation  coefficient 

A/D 

Analog-to-digital 

FIR 

Finite  impulse  response 

PSF 

Poisson  summation  formula 

TC 

Two’s-complement 

TCC 

Two’s-complement  chopping 

TCR 

Two’s-complement  roundoff 

ABBREVIATIONS 

l.s.b. 

Least  significant  bit 

p.d.f. 

Probability  density  function 

r.v. 

Random  variable 

SYMBOLS 

Equation 

C 

Auto-  and  cross-covariances 

(8) 

E 

Denotes  expectation  or  expected  value 

(5) 

I 

Indicator  function 

(27) 

K 

Multiplier  input  word  length  in  bits 

L 

Multiplication  coefficient  word  length  in  bits 

M 

Multiplier  output  word  length  in  bits  (result  of  roundoff  or  chopping) 

N 

Number  of  bits  discarded  in  the  roundoff  or  chopping  operation 
(K,L,M  do  not  include  sign  bit.  See  Fig.  2.) 

P 

Probability 

(64) 

Q 

Characteristic  function  or  Fourier  Transform  of  a  continuous 
function 

(65) 

R 

Auto-  and  cross-correlation 

(9) 

S 

Power  spectral  density 

(12) 

T 

Time  interval  between  consecutive  data  samples 

(1) 
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u 

Number  of  FIR  filter  taps 

(4) 

a 

Multiplication  coefficient  value 

(2) 

h 

Desired  or  ideal  filter  output  sequence 

(I) 

e 

Integer  form  of  the  multiplication  error  e 

(58) 

i 

v-1 

(12) 

k 

Integer  form  of  the  filter  input  sequence  x 

(46) 

k 

Integer  form  of  the  coefficient  a 

(47) 

9 

Quantization  step  size 

(64) 

r 

Remainder 

(35) 

u 

Integer  form  of  the  desired  filter  output  sequence  b 

(31) 

w 

Additive  noise  contribution  to  filter  output  sequence  y 

(1) 

X 

Filter  or  multiplier  input  data  sequence 

(2) 

y 

Filter  o>  multiplier  output  data  sequence 

(1) 

X 

Analog  waveform  used  for  quantizer  input 

* 

Denotes  machine  representation  (e.g..  b*>:  also,  complex  conjugate 

(28) 

(A-9) 

A 

Change  in  a  statistic  imposed  by  the  multiplication  errors 

(6) 

6 

Number  of  consecutive  zeroes  in  the  l.s.b.  positions  in  the  TC 
representation  of  the  coefficient  value 

(51) 

€ 

Multiplication  error  sequence 

(3) 

X 

Coefficient  related  odd  integer 

(104) 

M 

Mean  value 

(4) 

V 

Effective  word  length  in  bits  of  the  multiplication  error 

154) 

P 

Correlation  coefficient  (Ip)  <  1.0) 

a 

Standard  deviation;  with  no  subscript  o  denotes  o  j 

(5) 

LO 

Radian  frequency 

(12) 

